Assessment of the nano-mechanical properties of healthy and atherosclerotic coronary arteries by atomic force microscopy.

Nano-indentation techniques might be better equipped to assess the heterogeneous material properties of plaques than macroscopic methods but there are no bespoke protocols for this kind of material testing for coronary arteries. Therefore, we developed a measurement protocol to extract mechanical properties from healthy and atherosclerotic coronary artery tissue sections. Young's modulus was derived from force-indentation data. Metrics of collagen fibre density were extracted from the same tissue, and the local material properties were co-registered to the local collagen microstructure with a robust framework. The locations of the indentation were retrospectively classified by histological category (healthy, plaque, lipid-rich, fibrous cap) according to Picrosirius Red stain and adjacent Hematoxylin & Eosin and Oil-Red-O stains. Plaque tissue was softer (p < 0.001) than the healthy coronary wall. Areas rich in collagen within the plaque (fibrous cap) were significantly (p < 0.001) stiffer than areas poor in collagen/lipid-rich, but less than half as stiff as the healthy coronary media. Young's moduli correlated (Pearson's ρ = 0.53, p < 0.05) with collagen content. Atomic force microscopy (AFM) is capable of detecting tissue stiffness changes related to collagen density in healthy and diseased cardiovascular tissue. Mechanical characterization of atherosclerotic plaques with nano-indentation techniques could refine constitutive models for computational modelling.

A physics-based machine learning technique rapidly reconstructs the wall-shear stress and pressure fields in coronary arteries.

With the global rise of cardiovascular disease including atherosclerosis, there is a high demand for accurate diagnostic tools that can be used during a short consultation. In view of pathology, abnormal blood flow patterns have been demonstrated to be strong predictors of atherosclerotic lesion incidence, location, progression, and rupture. Prediction of patient-specific blood flow patterns can hence enable fast clinical diagnosis. However, the current state of art for the technique is by employing 3D-imaging-based Computational Fluid Dynamics (CFD). The high computational cost renders these methods impractical. In this work, we present a novel method to expedite the reconstruction of 3D pressure and shear stress fields using a combination of a reduced-order CFD modelling technique together with non-linear regression tools from the Machine Learning (ML) paradigm. Specifically, we develop a proof-of-concept automated pipeline that uses randomised perturbations of an atherosclerotic pig coronary artery to produce a large dataset of unique mesh geometries with variable blood flow. A total of 1,407 geometries were generated from seven reference arteries and were used to simulate blood flow using the CFD solver Abaqus. This CFD dataset was then post-processed using the mesh-domain common-base Proper Orthogonal Decomposition (cPOD) method to obtain Eigen functions and principal coefficients, the latter of which is a product of the individual mesh flow solutions with the POD Eigenvectors. Being a data-reduction method, the POD enables the data to be represented using only the ten most significant modes, which captures cumulatively greater than 95% of variance of flow features due to mesh variations. Next, the node coordinate data of the meshes were embedded in a two-dimensional coordinate system using the t-distributed Stochastic Neighbor Embedding (t-SNE) algorithm. The reduced dataset for t-SNE coordinates and corresponding vector of POD coefficients were then used to train a Random Forest Regressor (RFR) model. The same methodology was applied to both the volumetric pressure solution and the wall shear stress. The predicted pattern of blood pressure, and shear stress in unseen arterial geometries were compared with the ground truth CFD solutions on "unseen" meshes. The new method was able to reliably reproduce the 3D coronary artery haemodynamics in less than 10 s.

Considering the Influence of Coronary Motion on Artery-Specific Biomechanics Using Fluid-Structure Interaction Simulation.

The endothelium in the coronary arteries is subject to wall shear stress and vessel wall strain, which influences the biology of the arterial wall. This study presents vessel-specific fluid-structure interaction (FSI) models of three coronary arteries, using directly measured experimental geometries and boundary conditions. FSI models are used to provide a more physiologically complete representation of vessel biomechanics, and have been extended to include coronary bending to investigate its effect on shear and strain. FSI both without- and with-bending resulted in significant changes in all computed shear stress metrics compared to CFD (p = 0.0001). Inclusion of bending within the FSI model produced highly significant changes in Time Averaged Wall Shear Stress (TAWSS) + 9.8% LAD, + 8.8% LCx, - 2.0% RCA; Oscillatory Shear Index (OSI) + 208% LAD, 0% LCx, + 2600% RCA; and transverse wall Shear Stress (tSS) + 180% LAD, + 150% LCx and + 200% RCA (all p < 0.0001). Vessel wall strain was homogenous in all directions without-bending but became highly anisotropic under bending. Changes in median cyclic strain magnitude were seen for all three vessels in every direction. Changes shown in the magnitude and distribution of shear stress and wall strain suggest that bending should be considered on a vessel-specific basis in analyses of coronary artery biomechanics.

JAG1-NOTCH4 mechanosensing drives atherosclerosis.

Endothelial cell (EC) sensing of disturbed blood flow triggers atherosclerosis, a disease of arteries that causes heart attack and stroke, through poorly defined mechanisms. The Notch pathway plays a central role in blood vessel growth and homeostasis, but its potential role in sensing of disturbed flow has not been previously studied. Here, we show using porcine and murine arteries and cultured human coronary artery EC that disturbed flow activates the JAG1-NOTCH4 signaling pathway. Light-sheet imaging revealed enrichment of JAG1 and NOTCH4 in EC of atherosclerotic plaques, and EC-specific genetic deletion of Jag1 (Jag1ECKO) demonstrated that Jag1 promotes atherosclerosis at sites of disturbed flow. Mechanistically, single-cell RNA sequencing in Jag1ECKO mice demonstrated that Jag1 suppresses subsets of ECs that proliferate and migrate. We conclude that JAG1-NOTCH4 sensing of disturbed flow enhances atherosclerosis susceptibility by regulating EC heterogeneity and that therapeutic targeting of this pathway may treat atherosclerosis.

Mechanosensitive pathways are regulated by mechanosensitive miRNA clusters in endothelial cells.

Shear stress is known to affect many processes in (patho-) physiology through a complex, multi-molecular mechanism, termed mechanotransduction. The sheer complexity of the process has raised questions how mechanotransduction is regulated. Here, we comprehensively evaluate the literature about the role of small non-coding miRNA in the regulation of mechanotransduction. Regulation of mRNA by miRNA is rather complex, depending not only on the concentration of mRNA to miRNA, but also on the amount of mRNA competing for a single mRNA. The only mechanism to counteract the latter factor is through overarching structures of miRNA. Indeed, two overarching structures are present miRNA families and miRNA clusters, and both will be discussed in details, regarding the latest literature and a previous conducted study focussed on mechanotransduction. Both the literature and our own data support a new hypothesis that miRNA-clusters predominantly regulate mechanotransduction, affecting 65% of signalling pathways. In conclusion, a new and important mode of regulation of mechanotransduction is proposed, based on miRNA clusters. This finding implicates new avenues for treatment of mechanotransduction and atherosclerosis.

Considerations for analysis of endothelial shear stress and strain in FSI models of atherosclerosis.

Atherosclerosis is a lipid driven chronic inflammatory disease that is characterized by the formation of plaques at predilection sites. These predilection sites (side branches, curved segments, and bifurcations) have often been associated with disturbed shear stress profiles. However, in addition to shear stress, endothelial cells also experience artery wall strain that could contribute to atherosclerosis progression. Herein, we describe a method to accurately obtain these shear stress and strain profiles. We developed a fluid-structure interaction (FSI) framework for modelling arteries within a commercially available package (Abaqus, version 6.14) that included known prestresses (circumferential, axial and pressure associated). In addition, we co-registered 3D histology to a micro-CT-derived 3D reconstruction of an atherosclerotic carotid artery from a cholesterol-fed ApoE-/- mouse to include the spatial distribution of lipids within a subject-specific model. The FSI model also incorporated a nonlinear hyperelastic material model with regionally-varying properties that distinguished between healthy vessel wall and plaque. FSI predicted a lower shear stress than CFD (~-12%), but further decreases in plaque regions with softer properties (~-24%) were dependent on the approach used to implement the prestresses in the artery wall. When implemented with our new hybrid approach (zero prestresses in regions of lipid deposition), there was significant heterogeneity in endothelial shear stress in the atherosclerotic artery due to variations in stiffness and, in turn, wall strain. In conclusion, when obtaining endothelial shear stress and strain in diseased arteries, a careful consideration of prestresses is necessary. This paper offers a way to implement them.

Targeted Molecular Iron Oxide Contrast Agents for Imaging Atherosclerotic Plaque.

Overview: Cardiovascular disease remains a leading cause of death worldwide, with vulnerable plaque rupture the underlying cause of many heart attacks and strokes. Much research is focused on identifying an imaging biomarker to differentiate stable and vulnerable plaque. Magnetic Resonance Imaging (MRI) is a non-ionising and non-invasive imaging modality with excellent soft tissue contrast. However, MRI has relatively low sensitivity (micromolar) for contrast agent detection compared to nuclear imaging techniques. There is also an increasing emphasis on developing MRI probes that are not based on gadolinium chelates because of increasing concerns over associated systemic toxicity and deposits1. To address the sensitivity and safety concerns of gadolinium this project focused on the development of a high relaxivity probe based on superparamagnetic iron oxide nanoparticles for the imaging of atherosclerotic plaque with MRI. With development, this may facilitate differentiating stable and vulnerable plaque in vivo.Aim: To develop a range of MRI contrast agents based on superparamagnetic iron oxide nanoparticles (SPIONs), and test them in a murine model of advanced atherosclerosis. Methods: Nanoparticles of four core sizes were synthesised by thermal decomposition and coated with poly(maleicanhydride-alt-1-octadecene) (PMAO), poly(ethyleneimine) (PEI) or alendronate, then characterised for core size, hydrodynamic size, surface potential and relaxivity. On the basis of these results, one candidate was selected for further studies. In vivo studies using 10 nm PMAO-coated SPIONs were performed in ApoE-/- mice fed a western diet and instrumented with a perivascular cuff on the left carotid artery. Control ApoE-/- mice were fed a normal chow diet and were not instrumented. Mice were scanned on a 3T MR scanner (Philips Achieva) with the novel SPION contrast agent, and an elastin-targeted gadolinium agent that was shown previously to enable visualisation of plaque burden. Histological analysis was undertaken to confirm imaging findings through staining for macrophages, CX3CL1, elastin, tropoelastin, and iron. Results: The lead SPION agent consisted of a 10 nm iron oxide core with poly(maleicanhydride-alt-1-octadecene), (-36.21 mV, r2 18.806 mmol-1/s-1). The irregular faceting of the iron oxide core resulted in high relaxivity and the PMAO provided a foundation for further functionalisation on surface -COOH groups. The properties of the contrast agent, including the negative surface charge and hydrodynamic size, were designed to maximise circulation time and evade rapid clearance through the renal system or phagocytosis. In vitro testing showed that the SPION agent was non-toxic. In vivo results show that the novel contrast agent accumulates in similar vascular regions to a gadolinium-based contrast agent (Gd-ESMA) targeted to elastin, which accumulates in plaque. There was a significant difference in SPION signal between the instrumented and the contralateral non-instrumented vessels in diseased mice (p = 0.0411, student's t-test), and between the instrumented diseased vessel and control vessels (p = 0.0043, 0.0022, student's t-test). There was no significant difference between the uptake of either contrast agent between stable and vulnerable plaques (p = 0.3225, student's t-test). Histological verification was used to identify plaques, and Berlin Blue staining confirmed the presence of nanoparticle deposits within vulnerable plaques and co-localisation with macrophages. Conclusion: This work presents a new MRI contrast agent for atherosclerosis which uses an under-explored surface ligand, demonstrating promising properties for in vivo behaviour, is still in circulation 24 hours post-injection with limited liver uptake, and shows good accumulation in a murine plaque model.

Shear Stress Estimated by Quantitative Coronary Angiography Predicts Plaques Prone to Progress and Cause Events.

OBJECTIVES: This study examined the value of endothelial shear stress (ESS) estimated in 3-dimensional quantitative coronary angiography (3D-QCA) models in detecting plaques that are likely to progress and cause events. BACKGROUND: Cumulative evidence has shown that plaque characteristics and ESS derived from intravascular ultrasound (IVUS)-based reconstructions enable prediction of lesions that will cause cardiovascular events. However, the prognostic value of ESS estimated by 3D-QCA in nonflow limiting lesions is yet unclear. METHODS: This study analyzed baseline virtual histology (VH)-IVUS and angiographic data from 28 lipid-rich lesions (i.e., fibroatheromas) that caused major adverse cardiovascular events or required revascularization (MACE-R) at 5-year follow-up and 119 lipid-rich plaques from a control group that remained quiescent. The segments studied by VH-IVUS at baseline were reconstructed using 3D-QCA software. In the obtained geometries, blood flow simulation was performed, and the pressure gradient across the lipid-rich plaque and the mean ESS values in 3-mm segments were estimated. The additive value of these hemodynamic indexes in predicting MACE-R beyond plaque characteristics was examined. RESULTS: MACE-R lesions were longer, had smaller minimum lumen area, increased plaque burden (PB), were exposed to higher ESS, and exhibited a higher pressure gradient. In multivariable analysis, PB (hazard ratio: 1.08; p = 0.004) and the maximum 3-mm ESS value (hazard ratio: 1.11; p = 0.001) were independent predictors of MACE-R. Lesions exposed to high ESS (>4.95 Pa) with a high-risk anatomy (minimal lumen area <4 mm2 and PB >70%) had a higher MACE-R rate (53.8%) than those with a low-risk anatomy exposed to high ESS (31.6%) or those exposed to low ESS who had high- (20.0%) or low-risk anatomy (7.1%; p < 0.001). CONCLUSIONS: In the present study, 3D-QCA-derived local hemodynamic variables provided useful prognostic information, and, in combination with lesion anatomy, enabled more accurate identification of MACE-R lesions.

Reliable in vivo intravascular imaging plaque characterization: A challenge unmet.

Intravascular imaging has enabled in vivo assessment of coronary artery pathology and detection of plaque characteristics that are associated with increased vulnerability. Prospective invasive imaging studies of coronary atherosclerosis have demonstrated that invasive imaging modalities can detect lesions that are likely to progress and cause cardiovascular events and provided unique insights about atherosclerotic evolution. However, despite the undoubted value of the existing imaging techniques in clinical and research arenas, all the available modalities have significant limitations in assessing plaque characteristics when compared with histology. Hybrid/multimodality intravascular imaging appears able to overcome some of the limitations of standalone imaging; however, there are only few histology studies that examined their performance in evaluating plaque pathobiology. In this article, we review the evidence about the efficacy of standalone and multi-modality/hybrid intravascular imaging in assessing plaque morphology against histology, highlight the advantages and limitations of the existing imaging techniques and discuss the future potential of emerging imaging modalities in the study of atherosclerosis.

Disturbed flow induces a sustained, stochastic NF-κB activation which may support intracranial aneurysm growth in vivo.

Intracranial aneurysms are associated with disturbed velocity patterns, and chronic inflammation, but the relevance for these findings are currently unknown. Here, we show that (disturbed) shear stress induced by vortices is a sufficient condition to activate the endothelial NF-kB pathway, possibly through a mechanism of mechanosensor de-activation. We provide evidence for this statement through in-vitro live cell imaging of NF-kB in HUVECs exposed to different flow conditions, stochastic modelling of flow induced NF-kB activation and induction of disturbed flow in mouse carotid arteries. Finally, CFD and immunofluorescence on human intracranial aneurysms showed a correlation similar to the mouse vessels, suggesting that disturbed shear stress may lead to sustained NF-kB activation thereby offering an explanation for the close association between disturbed flow and intracranial aneurysms.

Plaque burden influences accurate classification of fibrous cap atheroma by in vivo optical coherence tomography in a porcine model of advanced coronary atherosclerosis.

AIMS: In vivo validation of coronary optical coherence tomography (OCT) against histology and the effects of plaque burden (PB) on plaque classification remain unreported. We aimed to investigate this in a porcine model with human-like coronary atherosclerosis. METHODS AND RESULTS: Five female Yucatan D374Y-PCSK9 transgenic hypercholesterolaemic minipigs were implanted with a coronary shear-modifying stent to induce advanced atherosclerosis. OCT frames (n=201) were obtained 34 weeks after implantation. Coronary arteries were perfusion-fixed, serially sectioned and co-registered with OCT using a validated algorithm. Lesions were adjudicated using the Virmani classification and PB assessed from histology. OCT had a high sensitivity, but modest specificity (92.9% and 74.6%), for identifying fibrous cap atheroma (FCA). The reduced specificity for OCT was due to misclassification of plaques with histologically defined pathological intimal thickening (PIT) as FCA (46.1% of the frames with histological PIT were misclassified). PIT lesions misclassified as FCA by OCT had a statistically higher PB than in other OCT frames (median 32.0% versus 13.4%; p<0.0001). Misclassification of PIT lesions by OCT occurred when PB exceeded approximately 20%. CONCLUSIONS: Compared with histology, in vivo OCT classification of FCA had high sensitivity but reduced specificity due to misclassification of PITs with high PB.

Connexin40 controls endothelial activation by dampening NFκB activation.

Connexins are proteins forming gap junction channels for intercellular communication. Connexin40 (Cx40) is highly expressed by endothelial cells (ECs) of healthy arteries but this expression is lost in ECs overlying atherosclerotic plaques. Low/oscillatory shear stress observed in bends and bifurcations of arteries is atherogenic partly through activation of the pro-inflammatory NFκB pathway in ECs. In this study, we investigated the relation between shear stress, Cx40 and NFκB. Shear stress-modifying casts were placed around carotid arteries of mice expressing eGFP under the Cx40 promoter (Cx40+/eGFP ). We found that Cx40 expression is decreased in carotid regions of oscillatory shear stress but conserved in high and low laminar shear stress regions. These results were confirmed in vitro. Using phage display, we retrieved a binding motif for the intracellular regulatory Cx40 C-terminus (Cx40CT), i.e. HS[I, L, V][K, R]. One of the retrieved peptides (HSLRPEWRMPGP) showed a 58.3% homology with amino acids 5-to-16 of IκBα, a member of the protein complex inhibiting NFκB activation. Binding of IκBα (peptide) and Cx40 was confirmed by crosslinking and en face proximity ligation assay on carotid arteries. TNFα-induced nuclear translocation of NFκB in ECs was enhanced after reducing Cx40 with siRNA. Transfection of HeLa cells with either full-length Cx40 or Cx40CT demonstrated that Cx40CT was sufficient for inhibition of TNFα-induced NFκB phosphorylation. Finally, Tie2CreTgCx40fl/flApoe-/- mice showed exaggerated shear stress-induced atherosclerosis and enhanced NFκB nuclear translocation. Our data show a novel functional IκBα-Cx40 interaction that may be relevant for the control of NFκB activation by shear stress in atherogenesis.

Interferon Regulatory Factor 5 Controls Necrotic Core Formation in Atherosclerotic Lesions by Impairing Efferocytosis.

BACKGROUND: Myeloid cells are central to atherosclerotic lesion development and vulnerable plaque formation. Impaired ability of arterial phagocytes to uptake apoptotic cells (efferocytosis) promotes lesion growth and establishment of a necrotic core. The transcription factor interferon regulatory factor (IRF)-5 is an important modulator of myeloid function and programming. We sought to investigate whether IRF5 affects the formation and phenotype of atherosclerotic lesions. METHODS: We investigated the role of IRF5 in atherosclerosis in 2 complementary models. First, atherosclerotic lesion development in hyperlipidemic apolipoprotein E-deficient (ApoE-/-) mice and ApoE-/- mice with a genetic deletion of IRF5 (ApoE-/-Irf5-/-) was compared and then lesion development was assessed in a model of shear stress-modulated vulnerable plaque formation. RESULTS: Both lesion and necrotic core size were significantly reduced in ApoE-/-Irf5-/- mice compared with IRF5-competent ApoE-/- mice. Necrotic core size was also reduced in the model of shear stress-modulated vulnerable plaque formation. A significant loss of CD11c+ macrophages was evident in ApoE-/-Irf5-/- mice in the aorta, draining lymph nodes, and bone marrow cell cultures, indicating that IRF5 maintains CD11c+ macrophages in atherosclerosis. Moreover, we revealed that the CD11c gene is a direct target of IRF5 in macrophages. In the absence of IRF5, CD11c- macrophages displayed a significant increase in expression of the efferocytosis-regulating integrin-ß3 and its ligand milk fat globule-epidermal growth factor 8 protein and enhanced efferocytosis in vitro and in situ. CONCLUSIONS: IRF5 is detrimental in atherosclerosis by promoting the maintenance of proinflammatory CD11c+ macrophages within lesions and controlling the expansion of the necrotic core by impairing efferocytosis.

Zebrafish Model for Functional Screening of Flow-Responsive Genes.

OBJECTIVE: Atherosclerosis is initiated at branches and bends of arteries exposed to disturbed blood flow that generates low shear stress. This mechanical environment promotes lesions by inducing endothelial cell (EC) apoptosis and dysfunction via mechanisms that are incompletely understood. Although transcriptome-based studies have identified multiple shear-responsive genes, most of them have an unknown function. To address this, we investigated whether zebrafish embryos can be used for functional screening of mechanosensitive genes that regulate EC apoptosis in mammalian arteries. APPROACH AND RESULTS: First, we demonstrated that flow regulates EC apoptosis in developing zebrafish vasculature. Specifically, suppression of blood flow in zebrafish embryos (by targeting cardiac troponin) enhanced that rate of EC apoptosis (≈10%) compared with controls exposed to flow (≈1%). A panel of candidate regulators of apoptosis were identified by transcriptome profiling of ECs from high and low shear stress regions of the porcine aorta. Genes that displayed the greatest differential expression and possessed 1 to 2 zebrafish orthologues were screened for the regulation of apoptosis in zebrafish vasculature exposed to flow or no-flow conditions using a knockdown approach. A phenotypic change was observed in 4 genes; p53-related protein (PERP) and programmed cell death 2-like protein functioned as positive regulators of apoptosis, whereas angiopoietin-like 4 and cadherin 13 were negative regulators. The regulation of perp, cdh13, angptl4, and pdcd2l by shear stress and the effects of perp and cdh13 on EC apoptosis were confirmed by studies of cultured EC exposed to flow. CONCLUSIONS: We conclude that a zebrafish model of flow manipulation coupled to gene knockdown can be used for functional screening of mechanosensitive genes in vascular ECs, thus providing potential therapeutic targets to prevent or treat endothelial injury at atheroprone sites.

Disturbed Cyclical Stretch of Endothelial Cells Promotes Nuclear Expression of the Pro-Atherogenic Transcription Factor NF-κB.

Exposure of endothelial cells to low and multidirectional blood flow is known to promote a pro-atherogenic phenotype. The mechanics of the vessel wall is another important mechano-stimulus within the endothelial cell environment, but no study has examined whether changes in the magnitude and direction of cell stretch can be pro-atherogenic. Herein, we developed a custom cell stretching device to replicate the in vivo stretch environment of the endothelial cell and examined whether low and multidirectional stretch promote nuclear translocation of NF-κB. A fluid-structure interaction model of the device demonstrated a nearly uniform strain within the region of cell attachment and a negligible magnitude of shear stress due to cyclical stretching of the cells in media. Compared to normal cyclical stretch, a low magnitude of cyclical stretch or no stretch caused increased expression of nuclear NF-κB (p = 0.09 and p < 0.001, respectively). Multidirectional stretch also promoted significant nuclear NF-κB expression, comparable to the no stretch condition, which was statistically higher than the low (p < 0.001) and normal (p < 0.001) stretch conditions. This is the first study to show that stretch conditions analogous to atherogenic blood flow profiles can similarly promote a pro-atherogenic endothelial cell phenotype, which supports a role for disturbed vessel wall mechanics as a pathological cell stimulus in the development of advanced atherosclerotic plaques.

Influence of shear stress magnitude and direction on atherosclerotic plaque composition.

The precise flow characteristics that promote different atherosclerotic plaque types remain unclear. We previously developed a blood flow-modifying cuff for ApoE-/- mice that induces the development of advanced plaques with vulnerable and stable features upstream and downstream of the cuff, respectively. Herein, we sought to test the hypothesis that changes in flow magnitude promote formation of the upstream (vulnerable) plaque, whereas altered flow direction is important for development of the downstream (stable) plaque. We instrumented ApoE-/- mice (n = 7) with a cuff around the left carotid artery and imaged them with micro-CT (39.6 µm resolution) eight to nine weeks after cuff placement. Computational fluid dynamics was then performed to compute six metrics that describe different aspects of atherogenic flow in terms of wall shear stress magnitude and/or direction. In a subset of four imaged animals, we performed histology to confirm the presence of advanced plaques and measure plaque length in each segment. Relative to the control artery, the region upstream of the cuff exhibited changes in shear stress magnitude only (p < 0.05), whereas the region downstream of the cuff exhibited changes in shear stress magnitude and direction (p < 0.05). These data suggest that shear stress magnitude contributes to the formation of advanced plaques with a vulnerable phenotype, whereas variations in both magnitude and direction promote the formation of plaques with stable features.